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ABSTRACT 


A  study  was  conducted  of  the  amplitude  modulation  of  an 
acoustic  signal  forward  scattered  from  a  wind-generated, 
model  ocean  surface.   The  demodulated  spectrum  of  the  surface 
scattered  sound,  was  compared  with  the  model  ocean  spectrum 
while  varying  the  acoustic  frequency  and  angle  of  incidence. 
Several  methods  of  enhancing  signal  detection  through  a 
knowledge  of  the  ocean  spectrum  were  studied.   The  most 
promising  of  these  shows  potential  gains  of  the  order  of 
three  decibels  when  the  received  signal  is  weak  and 
fluctuating  greatly,  as  at  a  Lloyd  Mirror  minimum. 
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I.   INTRODUCTION 

It  has  long  been  known  that  received  signals  from 
distant  underwater  sound  sources  vary  in  time.   The  effect 
has  been  described  as  "builds  and  fades"  of  the  signal. 
Though  this  effect  is  automatically  interpreted  by  the 
human  ear,  which  listens  when  the  sound  is  present  and  tunes 
out  when  it  disappears  in  the  noise,  little  or  no  advantage 
of  the  phenomenon  has  been  taken  by  the  extensive  signal 
processing  equipment  available  today.   Though  nominally  far 
superior  to  the  ear  in  response-time  and  frequency  resolution, 
this  equipment  often  fails  to  achieve  even  parity  with  a 
human  observer  when  processing  a  time  varying,  or  modulated, 
signal.   This  is  due  primarily  to  the  electronic  equipment's 
narrow  bandwidth  and  the  fact  that  it  integrates  continuously 
rather  than  selectively. 

It  was  the  objective  of  this  research  to  investigate  some 
of  the  parameters  affecting  this  phenomenon  of  a  modulated 
acoustic  signal  with  particular  attention  to  determining,  for 
surface  affected  signals: 

1.  The  best  laboratory  modeling  techniques. 

2.  The  relation  of  the  modulation  frequency  to  the 
surface  wave  frequency. 

3.  The  variation  of  the  modulation  amplitude  and  spectrum 
with  frequency  and  source/receiver  directivity. 

4.  The  optimization  of  signal  processing  of  an  amplitude 
modulated  signal  by  digital  techniques. 
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In  short  an  attempt  was  made  to  discover  a  "temporal 
filtering"  scheme  which  could  be  used  in  concert  with 
existing  high  speed  frequency  filtering  techniques,  such  as 
the  Past  Fourier  transform,  to  enhance  signal  processing. 
The  concept  has  been  dubbed  "EASE",  Environmentally  Adjusted 
Signal  Enhancement. 


12 


II.   THEORY 

The  subject  of  scattering  of  sound  from  the  ocean 
surface  has  been  treated  by  several  authors  utilizing  nearly 
as  many  methods.   The  following  is  a  summary  of  Refs.l  and  2 
which  point  up  those  features  of  the  scattered  spectrum 
which  are  germane  to  understanding  the  relationships  between 
sound  fluctuations  in  the  upper  ocean  and  the  ocean  wave 
spectrum. 

As  Parkins  summarizes  the  phenomenon:   "When  a  mono- 
chromatic acoustic  wave  is  incident  on  the  surface  of  the 
ocean  the  envelope  of  reradiation  varies  in  time  or  fades." 

Assuming  that  the  ocean  is  generated  by  a  stationary 
homogeneous,  Gaussian  process,  and  that  it  is  smooth  enough 
so  that  multiple  scattering  can  be  neglected,  the  incident 
acoustic  wave  of  harmonic  time  variation  e  w  is  given  by: 


.+  -ki'r  iwt 
$.  -  $  e     e 


with 


k,  =  k,  x  +  k,  y  +  k.,  z 

1    lx     lyJ  lz 


r\  /S.  /\  /\ 

=  -y  [^x  +  my  +  nz] 


,+ 


0   is  an  amplitude  term,  I,    m,  and  n  are  the  direction 
cosines  and  A  is  the  wave  length  of  the  incident  acoustic 
wave.   Other  factors  are  defined  in  Figure  1. 
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1  =  s(x,y,t) 


ACOUSTIC 
ENERGY 


OBS.  PT. 


FIGURE  1.   Coordinate  System  for  Theoretical  Development 

Following  the  development  of  Beckmann  and  Spizzichino 
[Ref.  2]  which  stipulates  a  source  free,  isovelocity  ocean, 
and  using  the  Helmholtz  Integral  development  of  Stratton 
[Ref.  3] j  the  farfield  reradiated  pressure  may  be  expressed 
as : 


1H 


i(ut  -kR„)   *  +  T 

*r  =  2iR- K—    "   e    L  dxdy 

o         z      L 


with; 


=  K  x  +  K„y  +  K  z  =  k  -k\ 
x     y  z  o    1 


and  other  factors  again  as  defined  in  Figure  1. 

The  irradiated  area  has  been  assumed  to  be  an  LxL 
aperture  which  is  several  wavelengths  across. 

In  this  research  it  was  the  spectral  density  of  the 
reradiated  pressure  that  was  of  interest.   This  may  be 
obtained  by  first  forming  the  autocovariance  function, 


R(t)  =  <$r(x,y,t)  $r*(x,y,t+x)> 


The  autocovariance  function  having  been  found,  the  Wiener- 
Khlnchin  Theorem  states  : 

oo 

S(n)  =JT{R(t)}  =   /  R(t)  elfiT  dx 

where  Q   Is  the  frequency  of  the  reradiated  (or  scattered) 
sound. 

Continuing  with  Parkins,  and  assuming  that  the  surface 
irregularities  are  much  smaller  than  the  aperture,  the 
expression  for  the  spectral  density  of  the  scattered  sound 


15 


Is  as  given  below.   It  should  be  noted  at  this  point  that 
the  requirement  for  small  surface  irregularities  corresponds 
to  a  roughness  parameter, 


•yfg   =  K  o   =  — t—  [cos  6,  +  cos  9_] 


as  defined  by  Beckmann  and  Spizzichino,  much  less  than  one. 
a   is  the  ocean  RMS  wave  height,  X  the  acoustic  wave  length 
and  8,  and  9p  the  angles  of  incidence  and  reflection, 
respectively  (see  Section  IV. A. 4  for  a  further  discussion 

of ->/£")• 

The  expression  for  the  spectral  density  of  the  scattered 
sound  is : 

-K  2c2        16  sin2K  I   sin2K  m 
S(fi)  =  Be   z    [6(fi-u)  p^3 2L_ 

VV 


+  2£2K  2a2{[A»(K  ,K  v)]26[n-to-(K  2+K2)k(G)hl 


+    [A'(_K  -K  )]26[Q-a)+(Kv2+Kv2)!i(G)^]}] 


with:   G  the  acceleration  of  gravity 

A'  the  Neumann-Pierson  directional  wave  spectrum 
B   the  amplitude  term 
6   the  Dirac  delta 
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Parkins  summarizes  this  equation:   "....(this  equation) 
shows  the  signal  reradiated  by  a  wind  generated  sea  to  have 
a  specularly  reflected  component  occurring  at  the  incident 
frequency  to  and  two  scattered  components  which  are  Doppler 
shifted  by  equal  amounts  from  to." 

The  spectrum  of  the  reradiated  sound  is  thus  seen  to  be 
frequency  modulated  by  the  Doppler  shift  term  and  amplitude 
modulated  by  the  wave  spectrum.   It  is  the  amplitude 
modulation  by  the  ocean  waves  that  is  considered  in  this 
research. 
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Ill .   EXPERIMENTAL  SETUP 

This  research  was  carried  out  using  the  Naval  Postgraduate 
School  OCEAN  ACOUSTIC  WAVE  FACILITY  (OAWF)  with  data 
collection  and  analysis  accomplished  with  a  hybrid  computer 
setup  designed  by  the  Special  Projects  Group  of  NADC  and 
affectionately  referred  to  as  OPHELEA  (OCEAN  PHYSICS 
ENVIRONMENT  EFFECTS  ANALYZER).   All  other  equipment  utilized 
was  standard  commercial  electronic  gear  and  the  author  feels 
that  a  detailed  explanation  of  its  functions  and  capabilities 
is  unnecessary  (see  Equipment  List  -  Section  III.D). 

A .   OAWF 

This  facility  consists  of  the  unique  combination  of  a 
wind-water  tunnel  and  an  anechoic  tank  (see  Figures  2  and  3)- 
Up  to  five  3/4  HP  centrifugal  fans  (three  in  this  case) 
comprise  a  Sea  Exciter  (Figure  4)  which  forces  air  down  the 
wind-water  tunnel,  which  is  about  17  meters  in  length, 
generating  waves  which  empty  into  an  anechoic  tank.   The  air 
space  in  the  tunnel  was  standardized  at  14.25  cm  and  the 
cross  section  of  the  entire  tunnel  measures  about  1.2  meters 
square.   Observation  of  the  waves  at  various  points  along 
the  tunnel  is  permitted  by  glass  observation  windows  shown 
in  Figure  5.   Waves  entering  the  3x3x2.9  meter  tank  portion 
of  the  facility  are  constrained  from  diffracting  by  a 
plexiglass  restraining  barrier  and  are  terminated  by  a 
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Figure  4.   OAWF  Sea  Exciter 


Figure  5 . 

OAWF  Observation  Window 
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"beach"  to  prevent  reflections.   The  beach  spans  the  width 
of  the  tank  experimental  area  and  consists  of  a  screen 
covered  wedge  filled  with  stainless  steel  millings  which 
cause  incident  waves  to  dissipate  rather  than  reflect. 
The  tank  itself  is  lined  with  redwood  4x4' s  on  all  sides 
and  the  bottom  which  ensures  excellent  sound  absorption  at 
the  boundaries  as  pointed  out  in  Ref.  4.   These  redwood 
beams  are  oriented  with  a  corner  facing  toward  the  tank 
interior,  except  on  the  wall  bordering  the  waves  which  is 
flat,  to  further  increase  absorption.   A  movable  race  permits 
two-dimensional  positioning  of  transducers  to  within  one 
centimeter.   This  facility  models  the  real  ocean  on  a  scale 
of  about  50-1  and  is  shown  in  operation  in  Figures  6  through 
8. 

B.   OPHELEA 

The  OPHELEA  system  consists  of  essentially  three 
components  which  are  interfaced  to  provide  high  speed  analog 
to  digital  data  collection  and  digital  data  processing  with 
somewhat  slower  display.   The  design  was  developed  by  the 
Special  Projects  Section  of  the  Naval  Air  Development  Center 
in  conjunction  with  Pinkerton  Computer  Consultants  Inc. , 
of  Warminster,  Pennsylvania.   The  three  components  are  shown 
in  Figure  9  and  described  below. 

1.   Interdata  Model  70  Computer 

This  minicomputer  is  modularly  constructed  and 
features  LSI  and  other  state-of-the-art  computer  technology. 
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Figure  9.   The  OPHELEA  System  with  Computer  and  A/D  Converter 
in  Large  Cabinet  to  the  Left 
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It  Is  a  16  Bit  half-word  oriented  digital  computer  that  is 
FORTRAN  programmable  and  has  sufficient  memory  (32  thousand 
8-bit  bytes)  capacity  to  process  in  excess  of  8000  data 
points  with  an  average  size  program.   In  this  research  It 
was  programmed  to  perform  a  Fast  Fourier  Transform  on  input 
data  and  output  spectral  density  vs.  frequency.   Some 
statistical  analysis  was  also  performed.   This  minicomputer 
occupies  the  upper  half  of  the  large  cabinet  shown  in 
Figure  9. 

2.   Phoenix  Analog  to  Digital  Converters  Model  ADC  712 
Two  ADC  712  converters  may  be  used  separately  or 
simultaneously.   Each  A/D  converter  is  a  high  speed  device 
capable  of  encoding  up  to  ±10  volt  input  signals  into 
12  bits  (11  bits  +  sign)  of  data,  providing  resolution  to 
one  part  in  4095  (accuracy  of  5mv  in  20v).   The  maximum  rate 
of  conversion  is  6.5  microseconds  and  it  has  performed 
reliably  at  a  sampling  rate  of  up  to  300kHz.   The  A/D  is 
driven  by  a  "command  to  convert"  signal  consisting  of  a 
5  volt  positive  pulse  of  at  least  200  microseconds  duration 
which  is  supplied  by  an  external  oscillator.   Measurement 
is  accomplished  by  successive  comparisons  against  an  internal 
reference  voltage  after  which  the  value  is  converted  to  a 
digital  number.   The  A/D  converters  are  located  in  part  of 
the  lower  half  of  the  OPHELEA  cabinet. 
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3.   Texas  Instruments  Silent  700  Electronic 
Data  Terminal  Model  733 

The  TI  733  consists  of  a  keyboard  which  is  used  as 

a  programming  and  input/output  control  device,  a  printer 

which  outputs  data  or  answers  queries  of  the  computer  on 

heat  sensitive  paper  and  a  play  back  and  record  section 

which  is  used  to  read  in  computer  programs  and  record/ 

playback  data.   In  this  research,  due  to  storage  limitations, 

data  were  sometimes  taken  A/D,  read  out  to  tape  cassettes, 

then  reread  in  for  processing  after  the  appropriate  program 

had  been  loaded.   A  capability  of  interfacing  with  the 

IBM-360/67  CP/CMS  system,  allowing  direct  access  to  the  NPS 

computer  facility,  was  also  utilized  in  processing  some 

statistical  wave  height  information. 

C.   SEQUENCE  OP  OPERATION 

The  sequence  of  operation  used  most  often  in  this 
research  was  as  follows: 

1.  An  analog  signal  (varying  sound  pressure,  instan- 
taneous wave  height,  etc.)  was  input  to  the  A/D  converter 
which,  being  triggered  at  twice  the  Nyquist  rate  (a  value 
found  to  result  in  minimal  aliasing  effects),  converted  it 
to  digital  data. 

2.  Those  data  were  then  either: 

a.  directly  processed  by  the  appropriate  program 
which  had  been  installed  in  the  computer,  or, 

b.  read  out  to  tape  cassettes  for  later  analysis 
as  described  above. 
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3.   Reduced  data  and  desired  results  were  displayed  on 
the  paper  readout  of  the  terminal.   Figure  10  is  a  pictorial 
description  of  this  process. 

D.   EQUIPMENT  LIST 

In  order  to  avoid  needless  repetition  of  the  lengthy 
manufacturers'  names  for  electronic  equipment  throughout 
this  work,  a  number  of  abbreviations  will  be  used  instead. 
These  short  titles  will  be  employed  for  the  most  part  in 
block  diagrams  and  schematics  but  may  also  appear  in  the 
text.  The  abbreviations,  with  the  corresponding  full 
equipment  title/description  are  as  listed  below: 


ABBREVIATION 
CORRELATION  COMPUTER 

DC  POWER  SUPPLY 
G.R.  OSCILLATOR 

HP  PRECISION  AMP 

HP  POWER  AMP 

KH  FILTER- 3.3^2 
KH  FILTER-3550 
0' SCOPE 


MANUFACTURERS  TITLE 

Princeton  Applied  Research 
Model  101  Correlation  Function 
Computer 

Hewlett-Packard  Model  721A  DC 
Power  Supply 

General  Radio  Model  1312  Decade 
Oscillator 

Hewlett-Packard  Model  463A 

Precision  Amplifier 

Hewlett-Packard  Model  *167A 
Power  Amplifier 

Khron-Hite  Model  33^2  Filter 

Khron-Hite  Model  3350  Filter 

Tektronix,  Inc.  Model  5^5B 
Oscilloscope 
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ABBREVIATION 


PAR  PREAMP 


WAVETEK  142 
WAVETEK  144 
X-Y  RECORDER 


MANUFACTURERS  TITLE 

Prfeeeton  Applied  Research  Model  113 
PreAmpllfler  (Note:   HFRO  and  LFRO 
refer  to  the  High/Low  frequency  roll 
off  feature  of  the  preamp  —  a 
filtering  option) 

Wavetek  Model  142  HF  VCG  Generator 

Wavetek  Model  144  HF  VCG  Generator 

Varian  Aerograph  Model  F-100  Strip 
Chart  Recorder 
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Signal 
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Figure  10.   Data  Processing  Sequence 
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IV.   EXPERIMENTAL  PRQGEDUBEJ  AND  RESULTS 

A.   OAWP  CHARACTERISTICS 

In  order  both  to  ensure  the  validity  of  the  research  to 
follow  and  to  gain  an  insight  into  the  significant  parameters 
affecting  the  scattering  of  sound  from  a  real  or  model  ocean, 
several  preliminary  experiments  were  conducted.   In  the 
sections  that  follow  each  will  be  fully  discussed  with  setup 
diagrams,  procedural  explanations  and  results  included. 

1.   Frequency  Spectrum 

In  determining  the  frequency  spectrum  of  the  ocean 
wave  heights,  a  highly  linear  probe  designed  at  the  Civil 
Engineering  Department  of  Stanford  University  and  pictured 
in  Figures  11  and  12  was  utilized.   The  probe  is  capacitive 
in  nature,  with  a  covered  wire  as  the  sensor.   The  wire 
serves  as  one  electrode,  the  water  as  the  other  and  the 
insulation  is  the  dielectric  of  the  capacitive  element. 
The  probe  is  excited  by  a  2.4  kHz  sinusoidal  carrier 
(amplitude  4.5  volts  peak  to  peak).   When  the  probe  is 
immersed  in  the  water  wave  system  the  output  across  a 
capacitive  bridge  (Figure  13)  is  amplitude  modulated  by 
the  passing  waves  rising  and  falling  around  the  wire  and 
thereby  changing  the  effective  length  of  the  capacitor. 
This  AM  signal  is  demodulated  by  the  detector  shown  in 
Figure  14  and  the  resultant  time  varying  DC  signal,  a  replica 
of  the  instantaneous  wave  height  at  one  point,  is  frequency 
analyzed  by  the  OPHELEA  system. 
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Figure  11.   Stanford 

Wave  Probe 


Figure  12. 
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Figure    13.      Wave  Height   Capacitive  Bridge 
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As  noted  above  frequency  analysis  was  accomplished 
by  the  OPHELEA  system,  using  a  program  called  PING  1C ,  a 
copy  of  which  is  included  in  Appendix  A.   In  this  case  the 
program  called  for  taking  256  A/D  samples  at  a  sampling  rate 
of  20  Hz  and  computing  a  128-point  spectral  density  versus 
frequency  output.   These  parameters  provide  a  frequency 
resolution  of  greater  than  .08  Hz. 

When  taking  A/D  samples  at  sampling  frequencies  in 
excess  of  20  kHz,  a  sinusoidal  oscillator  which  is  half-wave 
rectified  by  computer  circuitry  was  used  to  provide  the 
"command-to- con vert"  pulse  train.   At  low  frequencies  however 
the  sinusoidal  pulse  duration  is  too  great  and  erroneous 
multiple  triggering  of  the  A/D  converter  results.   Therefore 
the  WAVETEK  144,  with  a  positive  pulse  waveform  and  with  its 
asymmetry  control  set  for  minimum  pulse  duration,  was  employed 
for  the  low  sampling  frequencies  required  for  finding  spectra 
of  wave  and  demodulated  scattered  sound.   A  block  diagram 
of  the  experimental  setup  is  shown  in  Figure  15  and  the 
normalized  spectral  density  plot  is  shown  in  Figure  16. 
This  plot  is  the  average  of  forty  runs  conducted  as  outlined 
above  (i.e.,  40-128  point  FFT  analyses). 

As  can  be  seen  from  the  figure,  the  predominant  OAWF 
ocean  frequency  is  about  2.5Hz  as  compared  with  a  typical 
actual  ocean  frequency  on  the  order  of  0.1Hz.   Though  neither 
ocean  is  monochromatic,  the  OAWF  spectrum  is  fairly  narrow. 
Its  bandwidth  is  about  0.6Hz  at  the  half  power  points  and 
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Surface   Wave   Spectrum 
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the  spectral  density  tends  to  fall  off  less  rapidly  toward 
higher  frequencies  than  does  the  actual  ocean  spectrum 
[Ref.  5].   This  is  due  to  the  greater  contribution  of 
capillary  waves  in  the  OAWF  waves.   The  irregularity  at 
about  2.4Hz  was  felt  at  first  to  be  merely  a  coincidental 
artifact  which  would  disappear  with  additional  averaging. 
However,  though  it  smoothed  some,  it  has  persisted.   The 
source  of  the  irregularity  may  have  been  lack  of  adequate 
fetch  for  wave  growth  to  an  equilibrium  spectrum.   It  is 
noted  that  a  similar  effect  was  observed  in  the  work  of 
Fowler  and  Scheibel  [Ref.  6]  at  short  fetch. 
2.   Probability  Density  Function 

As  Kinsman  has  noted  [Ref.  7]  the  probability  density 
function  associated  with  wind  driven  ocean  wave  heights  is 
nearly  Gaussian.   In  order  to  verify  that  the  OAWF  sea  met 
this  criterion,  wave  height  samples  were  collected  using  the 
setup  shown  in  Figure  15  of  the  previous  section  and  analysis 
was  accomplished  using  the  standard  IBM-36O/67  Library 
program  "HISTF" .   This  process  was  accomplished  as  follows: 

a.  Eight  groups  of  1024  wave  height  samples  were 
taken  from  the  previously  described  wave  probe  using  the 
OPHELEA  system  and  a  program  modestly  entitled  PERK  1A  and 
were  stored  in  the  minicomputer  memory. 

b.  The  data  were  then  read  out  to  tape  cassettes  in 
a  revised  format  compatible  with  the  IBM-36O  CP/CMS  system. 

c.  Data  were  converted  to  punch  cards  using  the 
CP/CMS  system  and  analyzed  using  "HISTF". 
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The  "HISTF"  program  outputs  various  statistical 
parameters  and  plots  a  histogram  of  the  input  data.   Each  of 
the  eight  data  groups  was  analyzed  and  the  results  averaged. 
A  graph  of  the  resulting  average  histogram,  as  well  as  a 
comparable  normalized  Gaussian  PDF,  is  shown  in  Figure  17. 
As  can  be  seen,  the  histogram  distribution  is  slightly  skewed 
and  displays  negative  kurtosis.   Though  not  exactly  agreeing 
with  any  one  of  Kinsman's  seas,  the  variations  are  slight  as 
is  shown  by  the  following  comparison: 

Kinsman  -  record  072  OAWF  -  3  Fans 

Skewness      +.092  -.0245 

Kurtosis      -.031  -. 35^9 

The  mean  is  also  seen  to  be  non-zero  but  is  less 
than  0.1  centimeter. 

3.   Autocorrelation  Function 

The  autocorrelation  function  of  the  ocean  wave 
heights  is  characteristically  that  of  band  limited  noise  and 
takes  the  shape  of  an  attenuated  cosine  [Ref.  8].   In  finding 
the  correlation  function  of  the  OAWF  three  fan  ocean, 
essentially  the  same  setup  was  used  as  in  the  previous  two 
sections.   However  the  output  from  the  wave  probe, 
characterizing  the  surface,  was  terminated,  not  in  the 
OPHELEA  system  as  before,  but  in  a  PAR  CORRELATOR.   This 
device  is  capable  of  computing  either  auto-  or  cross- 
correlation  functions  and  displaying  the  result  on  an 
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oscilloscope  or  a  recorder.   In  this  case  an  X-Y  RECORDER 
was  used.   A  diagram  of  the  experimental  setup  is  shown  in 
Figure  18. 

An  average  surface  wave  height  correlation  function 
for  time  lags,  t  =  0.0  to  2.0  seconds,  was  arrived  at  by 
superimposing  the  eleven  runs  printed  on  the  X-Y  RECORDER 
and  shown  in  Figure  19.   The  resultant  average  is  shown  in 
Figure  20  and  can  be  seen  to  follow  the  expected  form.   The 
rate  of  "attenuation"  is  initially  somewhat  less,  than 
measured  by  Kinsman,  with  the  OAWF  correlation  at  the  first 
peak  down  to  0.7  compared  to  0.50  for  Kinsman's  ocean  data. 
The  second  peak  correlation  conversely  has  dropped  to  about 
0.3  while  Kinsman's  C(x)  remains  at  about  0.4  after  time 
lags  of  two  ocean  periods. 

The  correlation  function  loses  its  consistency  after 
two  ocean  periods  as  may  be  seen  from  Figures  19  and  21. 
Figure  21  shows  the  function  computed  with  a  precomputation 
period  =  T  =  2.0  seconds  set  into  the  PAR  CORRELATOR.   With 
this  setting  the  correlation  function  is  computed  for 
x=2.0  to  4.0  seconds  or  in  effect  the  function  of  Figure  19 
is  continued  for  another  two  seconds.   For  this  longer  time 
lag,  it  is  difficult  to  see  any  repeatability  in  successive 
correlation  functions. 

While  the  OAWF  values  are  somewhat  different  from  those 
of  the  Kinsman  ocean,  the  similarity  in  shape  is  obvious. 
This  validation  of  the  ocean  model  will  be  of  importance 
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later  in  this  work  when  comparisons  are  made  with  the 
modulation  envelope  of  scattered  sound. 
4.   The  Roughness  Parameter 

The  previously  mentioned  roughness  parameter ,  -\fg    , 
is  defined  by  Beckmann  and  Spizzichino  [Ref.  2]  as: 


2TTCT 


-/g"  =  — —  [cos  61  +  cos  02] 


with      a   =  ocean  RMS  wave  height 

X   =  acoustic  wave  length 

6, ,02  =  angles  of  incidence  and  reflection 5  respectively 
(measured  from  the  normal  to  the  ocean  surface). 

It  is  easily  seen  that  a  "smooth"  condition  (-rfg~  <<   1)  can 
be  achieved  through  use  of  low  acoustic  frequencies  or  large 
angles  of  incidence  (shallow  grazing  angles)  as  well  as  with 
the  obvious  small  RMS  surface  wave  heights.   The  converse  is 
of  course  also  true  for  "rough"  conditions  (ijg   >>   1). 
yi"  was  the  factor  that  was  used  to  model  the  OAWF  ocean  to 
the  actual  ocean  in  later  research  and  an  example  of  typical 
sets  of  conditions  in  both  domains  yielding  the  same 
-\lg   =  .076  is  given  bel< 
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Wind 

RMS  wave  height 
Acoustic  Frequency 
Source/Rcvr  Depth 
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As  noted  above  -y/g"  is  the  same  for  both  sets  of 
parameters  and  thus  the  typical  ocean  conditions  noted  in 
the  right  column  can  be  reproduced  if  the  values  in  the  left 
column  are  employed  in  the  OAWF. 

B.   COMPARISON  OF  WAVE  HEIGHT  AND  DEMODULATED  SCATTERED 
SOUND  SPECTRA  AND  CORRELATION  FUNCTIONS 

As  noted  in  the  Theory  section  the  frequency  spectrum  of 
the  amplitude  modulation  of  the  scattered  sound  is  directly 
related  to  the  frequency  spectrum  of  the  wind  driven  waves. 
Hence  it  is  important  to  compare  these  frequency  spectra. 
Further  it  is  of  interest  to  compare  autocorrelation  functions 
in  order  to  gain  an  insight  into  signal  processing  possibili- 
ties.  In  addition  to  these  two  comparisons  the  cross 
correlation  function  was  computed. 

1.   Frequency  Spectra  Comparison 

Though  some  trials  described  later  in  this  work  were 
conducted  with  a  directional  source,  the  principal  study 
used  an  omnidirectional  source  and  receiver,  which  were  felt 
to  best  model  the  submarine/sonobuoy  scenario.   The  source 
was  a  locally-constructed  transducer  made  from  a  2"  (o.d.) 
Glennite  Ceramic  sphere  which  was  covered  with  several  coats 
of  neopreme.   This  source  transducer  will  be  referred  to 
hereafter  as  "OMNI".   The  receiver  was  a  standard  Atlantic 
Research  LC-10  hydrophone  and  its  response  was  immediately 
amplified  by  30  dB,  NUS  Corporation  FET  preamp  powered  by 
a  12  volt  DC  POWER  SUPPLY.   Both  the  LC-10/NUS  combination 
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and  the  OMNI,  as  well  as  a  directional  piston  type  transducer 
(F-27)  used  later,  are  shown  in  Figure  22. 

As  the  spectra  of  the  OAWF  sea  waves  had  already 
been  computed  from  an  average  of  forty  runs,  only  the 
average  sound  amplitude  modulation  spectrum  remained  to  be 
determined.   To  accomplish  this  the  scattered  AM  sound  was 
demodulated  by  the  detector  shown  in  Figure  14  and  the 
resulting  time-varying  DC  signal  sent  to  the  OPHELEA  system. 
A  photograph  of  the  wave  forms  of  the  AM  scattered  sound, 
the  demodulated  scattered  sound  and  the  instantaneous  wave 
height  are  shown  in  Figures  23  and  24.   A  complete  schematic 
of  the  experimental  setup  is  shown  in  Figure  25. 

Forty  runs  were  conducted  with  256  A/D  samples  of 
the  demodulated  acoustic  20  kHz  signal  being  taken  for  each 
run  at  a  sampling  frequency  of  2  0  Hz,  again  using  program 
PING  1C.   As  before,  a  128-point  spectral  density  versus 
frequency  readout  was  calculated  for  each  run  and  the 
resulting  forty  spectra  averaged  and  normalized  for  comparisons 
with  the  previously  determined  average  wave  spectrum.   These 
two  averaged,  normalized  spectra  are  plotted  in  Figure  26 
which  shows  the  obvious  similarity.   The  demodulated  sound 
spectrum  has  a  shape  that  is  remarkably  similar  to  its 
counterpart  though  its  peaks  are  about  0.2  Hz  less  than  the 
predominant  wave  frequencies.   This  difference  is  presumably 
due  to  the  surface  drift  velocity  [Ref.  9].   The  3  dB  bandwidth 
of  the  demodulated  sound  spectrum  is  about  0.5  Hz  as  compared 
with  0.6  Hz  for  the  wave  spectrum. 
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Figure  22. 


Hydrophone  (LC-10)  and  Transducers  Used 
in  this  Research 
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Figure  23. 

Amplitude  Modulated 
Sound  (top)  and 
Instantaneous  Wave 
Height 


Figure  2h  . 

Demodulated  AM  Sound  (top) 
and  Instantaneous  Wave 
Height 
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Figure    25.      Experimental   Setup    for  Determining  the   Spectrum 
of  Demodulated  Scattered  Sound 


51 


c 

o 
co 

■d 

^ 

■p 
-p 
cd 
o 
co 

■a 

cd 

o 

■a 

-p 
cd 

iH 

O 

s 

0) 
Q 


>> 

o 

o 

CD 

cd  -P 

3 

H    .i-, 

CT 

P  to 

$- 

Ll_ 

O  -H 

<D    (D 

CO 

>>  > 
O    Cd 

a; 

3  a) 

c  o 

0  cd 

73  CO 

cu 
•H  3= 

cd  O 

e 
o  c 

2  cd 


C\J 


bO 
•H 

I'. 


A1ISN3Q    IVdlOJJS  Q3ZI"lVWbON 


52 


These  frequency  spectra,  though  exhibiting  some  minor 
differences,  are  considered""tc --!_._.. -assent ially  the  same  and 
confirm  the  theoretically  expected  result. 
2 .   Autocorrelation  Function  Comparison 

The  same  procedure  was  used  in  determining  the 
demodulated  scattered  sound  autocorrelation  function  as  in 
the  case  of  the  wave  height  version.   The  output  of  the 
previous  setup  was  terminated  in  the  PAR  CORRELATOR  instead 
of  OPHELEA  and  recording  was  done  on  the  X-Y  RECORDER  as 
before.   Again  eleven  runs  were  conducted  for  delays 
t  =  0.0  to  2.0  seconds.   These  recordings  are  shown  in 
Figure  27.   An  average  autocorrelation  function  was  then 
fitted  to  this  record.   A  plot  of  both  records,  wave  height 
and  demodulated  scattered  sound,  is  shown  in  Figure  28. 
Though  the  ordinate  at  t  =  0  is  slightly  different  for  the 
two  functions  due  to  idiosyncracies  of  the  PAR  CORRELATOR 
and  X-Y  RECORDER,  it  is  apparent  that  the  results  are  nearly 
identical  for  the  two  cases.   Further,  it  is  evident  that 
significant  correlation  for  both  variables  exists  for  three 
"peaks",  that  is,  over  two  ocean  periods.   The  demodulated 
sound  correlation  function  also  becomes  confused  for  delays 
greater  than  2.0  seconds,  as  is  the  case  with  the  wave  height, 
as  seen  in  Figure  29.   It  should  be  noted  that  the  vertical 
scale  in  Figure  29  is  expanded  over  that  of  Figure  21  and  no 
inference  of  greater  correlation  should  be  drawn. 


53 


5H 


CD    CD 
>  -r- 

fO     CD 


x> 

CD 

•P 

cd 

rH 

3 

'O 

o 

6 

CD 

-P 

ax, 

M 

<M   -H 

O 

CO 

K 

CO 

C 

CD 

O 

> 

•H 

ctf 

-P^ 

o 

§ 

■a 

Pn 

a 

n 

•a 

o 

c 

•H 

? 

■P 

O 

cd  CO 

rH 

CO  'O 

Sn 

CO 

*H 

U 

O 

<D 

O  +3 

O  -P 

-p 

CTJ 

3 

o 

<  CO 

CO 

>a 

M  U 

cd 

cd 

?-< 

|5 

0 

m 

> 

o 

<  fe 

• 

CO 

C\J 

CO 

£-. 

3 

M 

•H 

Ua 

55 


->  -' 

I      CM 


1 

Isi  : 

cu 

1  s- 

1  CT> 

1 

o 

1 

56 


3.   Cross  Correlation  Function 

Using  the  two  setups ^described  in  the  autocorrelation 
section  above  the  two  wave  forms  were  cross  correlated, again 
using  the  PAR  CORRELATOR  and  X-Y  RECORDER.   The  same  sort  of 
averaging  process  was  then  performed  on  the  record  shown  in 
Figure  30  and  yielded  an  average  cross  correlation  function 
which  is  shown  in  Figure  31.   The  t  =  0  ordinate  varied 
considerably  between  positive  and  negative  values  of 
correlation  and  the  choice  of  a  near  zero  initial  value  was 
purely  arbitrary.   Again  three  peaks  are  seen  to  dominate, 
at  about  x  =  0.8,  1.25  and  1.7  seconds.   This  is  to  be 
expected  from  the  cross  correlation  function  of  two  signals 
which  both  exhibit  significant  correlation  over  two  ocean 
cycles  (or  have  three  "peaks"  in  their  autocorrelation 
functions) . 

C.   EFFECT  OF  FREQUENCY  AND  ANGLE  OF  INCIDENCE 
ON  SCATTERED  SOUND  MODULATION 

To  determine  whether  the  modulation  frequency  of  the 

scattered  sound  remained  essentially  constant  with  changing 

acoustic  frequency  and  varying  angle  of  incidence,  a  series 

of  trials  was  conducted  with  different  combinations  of  these 

two  parameters.   It  was  felt  at  this  stage  of  the  research 

that  there  might  be  some  advantage  in  using  a  directional 

hydrophone  to  make  the  detection. 
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1.   Changing  Frequency  and  Angle  of  Incidence 

In  these  tests  a  directional  transducer,  a  USRD 
Type  F-27,  was  mounted  in  a  trainable  bracket  which  permitted 
angles  of  incidence,  measured  from  the  surface  normal,  of 
from  0°  to  90°.   The  F-27  is  shown  in  Figure  22.   The 
scattered  sound  was  received  by  an  omnidirectional  LC-10 
hydrophone  which  was  positioned  at  the  same  depth  as  the 
F-27  and  2.75  meters  away.   To  achieve  different  angles  of 
incidence  while  retaining  specular  scatter  the  F-27  and 
LC-10  were  positioned  at  depths  of  0.99,  0.49  and  0.27  meters 
which  equate  to  angles  of  incidence  of  approximately  58°, 
70°  and  79°  respectively.   Prior  to  each  trial  the  F-27  was 
elevated  and  depressed,  with  the  surface  quiescent,  until 
the  maximum  response  ensured  that  the  observation  was  of 
specular  scatter.   Ten  runs,  using  program  PING  1C  (Appendix 
A),  were  conducted  at  each  position  with  256  A/D  samples 
being  taken  at  a  sampling  rate  of  20  Hz.   A  128-point  spectral 
density  versus  frequency  output  was  generated  as  in  the 
determination  of  previous  spectra.   Again  the  WAVETEK 
oscillator  provided  the  pulse  train  for  proper  A/D  converter 
triggering  at  low  frequencies.   A  diagram  of  this  setup  is 
shown  in  Figure  32  and  a  table  of  run  parameters  and 
particular  equipment  settings  is  found  in  Table  II. 

After  the  above  80  runs  were  completed,  the  ten 
spectra  from  each  angle/frequency  combination  were  averaged 
and  normalized  to  permit  comparison.   The  results  were  then 
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Figure  32.   Experimental  Setup  for  Determining  the  Effects 
of  Varying  Angle  of  Incidence  and  Frequency 
on  Scattered  Sound  Modulation 
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plotted  as  indicated  in  Table  I  below  with  results  compared 
for  constant  angle  of  incidence  and  variable  frequency  and 
vice  versa. 

TABLE  I 

■SUMMARY  OF  VARIABLE  ANGLE  OF 
INCIDENCE  AND  FREQUENCY  TRIALS 

Figure  Acoustic  Frequency   Angle  of  Incidence 

33  30,40  kHz  79° 

34  20,30,40  kHz  70° 

35  20,30,40  kHz  58° 

36  20  kHz  70°, 58° 

37  30  kHz  79°, 70°, 58° 

38  40  kHz  79°, 70°, 58° 

As  may  be  seen  from  Figures  33  through  38  the 
frequency  spectral  peaks  differ  by  less  than  0.3  Hz  in  all 
cases  except  Figure  35.   The  somewhat  wider  spread,  about 
0.5  Hz,  between  peaks  is  thought  to  be  due  to  greater 
contributions  from  the  side  lobes  at  that  angle  cf  incidence. 
Even  this  difference,  though  larger,  is  comparable  to  the 
3  dB  bandwidth  of  the  wave  frequency  spectrum. 

Based  on  this  work,  it  was  concluded  that  the 
modulation  spectra  of  the  scattered  sound  received  by  a 
point  hydrophone  are  similar  for  variations  in  acoustic 
frequency  or  angle  of  incidence.   It  is  felt  that  the 
differences  are  due  to  inadequate  averaging  and/or  the 
dependence  on  side  lobes. 
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2.   Modulation  Strength  with  Changing  Receiver  Angle 

The  object  of  this  set  of  trials  was  to  determine  if 
the  modulation  effect  could  be  increased  by  training  a 
directional  receiver  to  a  preferred  angle  of  reflection  with 
a  point  source.   Although  acoustic  reciprocity  would  be 
expected  to  hold,  in  this  case  control  would  be  at  the  option 
of  the  listener.   Specifically  it  was  felt  that  possibly  if 
the  receiver  hydrophone  were  directed  normal  to  the  modulating 
ocean  surface,  the  relative  modulation  of  the  received 
scattered  acoustic  signal  would  be  increased.   The  experi- 
mental setup  was  essentially  as  before;  however,  now  the 
P-27  was  the  receiver  and  the  previously  referred  to  OMNI  was 
the  source,  with  both  at  a  depth  of  1.0  meters.   The  three 
fan  OAWP  sea  was  again  the  environment  with  the  water  level 
standardized  as  before.   A  setup  diagram  Is  shown  in  Figure 
39.   Ten  runs  were  conducted  at  each  of  seven  angles  of 
incidence  ranging  from  direct  path,  through  specular  to 
normal  incidence.   The  acoustic  frequency  was  30  kHz  and 
amplifier  and  filter  settings  were  kept  constant  for  all 
runs  so  that  valid  comparisons  could  be  made.   Table  III 
shows  a  summary  of  these  trials. 

Each  of  the  seventy  runs  in  Table  III  consisted  of 
taking  256  A/D  samples  at  a  20  Hz  sampling  frequency  (using 
program  PINC  1C)  as  before;  then  averaging  the  spectral 
density  Information  for  each  angular  setting.   The  plot  of 
these  average  spectra  Is  presented  in  Figure  '10  which  shows 
log  spectral  density  plotted  versus  scattered  sound 
modulation  frequency. 

70 


TABLE  III 
SUMMARY  OP  VARIAB^^s^LE-JBEGEITSSR  TRIALS 

Run  Numbers  Angle  of  Incidence 

61-70  90°  (Direct) 

71-80  58°  (Specular) 

81-90  45° 

91-100  30° 

101-110  0°   (Normal 

111-120  25° 

121-130  75° 

Note:   Again  the  angle  of  Incidence  is  measured 
normal  to  the  surface. 

All  levels  have  been  plotted  without  adjustment 
so  that  the  amplitude  of  the  modulation  at  different  angles 
of  incidence  can  be  seen  simultaneously.   Summarizing: 

a.  The  greatest  modulation  strength  was  observed 
at  9  =  58°,  that  is  the  specular  direction. 

b.  Modulation  amplitude  declined  regularly  as  the 
angle  of  receiver  is  changed  in  either  direction  from  the 
specular. 

c.  Although  the  incoherent  component  (modulated 
amplitude)  may  possibly  be  increased  relative  to  the 
coherent  component,  the  modulation  amplitude  is  down  30  dB 
(re  specular)  with  the  receiver  trained  normal  to  the  sea 
surface  (0  =  0° ) . 

Knowing  that  the  amplitude  modulation  is  greatest  at 
or  near  the  specular  direction  might  provide  a  means  of 
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arriving  at  a  rough  fix  or  target  probability  area.   Assuming 
that  target  bearing  could  be  provided,  by  a  directional  hydro- 
phone and  that  the  local  sound  velocity  profile  is  known,  the 
approximate  range  could  be  found  by  applying  ray  path 
techniques  once  the  angle  of  greatest  modulation  (specular 
direction)  is  found.   This  scheme  would  of  course  require 
a  hydrophone  system  which  is  directional  in  both  azimuth 
and  elevation. 
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V.   ENHANCING  THE  TIME-VARYING  SIGNAL 

At  this  stage,  with  reasonable  assurance  that  the  OAWP 
sea  accurately  modeled  the  actual  ocean  and  that  the  modula- 
tion spectrum  of  the  scattered  sound  was  essentially  that  of 
the  surface  waves,  efforts  were  begun  toward  improving  signal 
processing  using  the  model  and  the  information  gained  to  date, 
If  the  ocean  wave  spectrum  were  of  a  single  fixed  frequency 
the  sound  would  fluctuate  at  that  same  single  fixed  frequency 
and  sound  signal  enhancement  would  then  be  easy.   It  would  be 
necessary  merely  to  identify  a  loud  point  in  the  received 
signal  and  then  to  listen  at  later  intervals  equal  to  the 
ocean  wave  period.   The  ocean  is  regrettably  not  this  simple 
as  is  obvious  from  the  respective  frequency  spectra  and 
therefore  alternative  or  modified  processing  means  had  to  be 
identified  and  implemented. 

A.   PROGRAMS  PERK  1A  AND  PERK  IB 

To  begin  to  explore  the  problem,  a  pair  of  computer 
programs  modestly  named  for  this  rookie  acoustician  were 
developed  which  would  accomplish  the  following: 

1.  PERK  1A  orders  the  computer  to  direct  the  A/D 
converter  to  take  a  group  of  samples  (say  64),  hereafter 
called  a  "BLOCK",  and  then  to  "pause". 

2.  The  "pause"  is  triggered  by  the  computer  taking  its 
last  ordered  sample  and  PERK  1A  then  shifts  the  computer 
into  a  Fortran  "Do  Loop"  which  has  an  Index  proportional  to 
the  desired  time  delay. 
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3.   When  the  "Do  Loop"  has  completed  its  preset  number 
of  operations  (and  the  delay  or  pause  has  been  completed)  the 
program  signals  the  A/D  to  take  another  BLOCK  of  samples. 
This  is  followed  by  another  "pause",  another  BLOCK,  etc. 

The  net  result  is  a  series  of  BLOCKS  of  A/D  samples 
taken  at  specified  intervals  along  the  modulated  signal 
which,  when  a  total  of  8192  data  points  have  been  taken  (the 
present  storage  limit  of  OPHELEA),  are  read  out  to  a  tape 
cassette.   This  data  tape  is  then  put  aside  as  the  second 
program,  further  modestly  entitled  PERK  IB,  is  loaded. 

PERK  IB  reads  the  time  series  data  of  the  previously 
recorded  sample  BLOCKS  back  into  memory  destructively  and 
performs  an  FFT  spectral  analysis  on  each  BLOCK.   The  program 
may  then  be  directed  to  output  all  of  the  complete  spectral 
analyses,  or  just  the  major  spectral  peak  in  each  analysis. 
Thus  a  sequence  of  successive  sound  spectral  density  peaks 
is  available,  each  separated  by  a  fixed  time  delay,  the 
"pause".   The  duration  of  the  "pause"  is  selected  to  provide 
the  desired  number  of  spectral  densities  of  predetermined 
frequency  resolution  during  an  "ocean  period".   An  "ocean 
period"  is  here  defined  as  the  reciprocal  of  the  peak 
frequency  in  the  ocean  spectrum.   The  effect  of  the  modulation 
caused  by  the  ocean  is  vividly  illustrated  by  Figure  4l,  which 
depicts  the  peak  sound  spectral  density  versus  time,  and  which 
shows  a  swing  of  about  18  dB  between  maxima  and  minima  of 
signal  strength.   Copies  of  PERK  1A  and  PERK  IB  may  be  found 
in  Appendix  A. 


76 


CD 
C\J 


O 

o 


o 

en 


o 

00 


o 

UJ 


ii 


o    go 

c_> 
o 


o 


O 
LT> 


CQ 


DO 


O 

o 


■  o 
«3- 


o 


O 

■H 
-P 
CQ 

O 

o 

< 

a> 

<u 
-p 
-p 
cd 
o 

C/J 

U 

cd 

O- 

cd 

o 

cu 
<& 

p 

•H 


w 
C 
o 

•H 
P 

cd 

p 
o 

H 


a  0) 
S  bO 
<;  U 
cd 
>>J 
P 

•H     I 

to 
C  H 

cy  cd 
P  C 
hO 
rH.H 
Cd  C/j 

p  H 
c_>  cd 

<D    (D 

cx:s 

CO  w 


i  o 

to 

CO 

C     • 

•H  tH 

hO  cd 

c  £ 

cd   hO 

L° 

JG  -h 

C\J 

O  00 

cu 

■H 


(jpaun)   AlISfOQ  !\AJlD3dS 


77 


In  this  work,  the  BLOCKS  almost  invariably  contained 
64  A/D  samples  and  the  "pause"  was  selected  to  provide  eight 
BLOCKS  per  ocean  period.   The  "following  is  a  sample  delay 
time  calculation  for  a  20  kHz  acoustic  signal,  modulated  by 
a  2.5  Hz  ocean  (ocean  period  0.4  sec.)  and  sampled  at  80  kHz. 

For  64  Samples/BLOCK  at  80  x  10^  samples/sec: 
sampling  time/BLOCK  =0.8  msec. 

For  8  BLOCKS/ocean  period:  sampling  time/ocean  period  = 
6.4  msec.  Therefore,  for  one  ocean  period,  total  delay  time 
=  400  msec  -  6.4  msec  =  393.6  msec. 

Delay  between  BLOCKS  =  393.6/8  =  49.2  msec. 

Total  data  collected  in  16  ocean  periods 

_  ,64  samples.,  ,8  BLOCKS    ■,  ,16  ocean  periods., 
1    BLOCK   nocean  periodM    run  ' 

8lQ?  samPles 
run 

For  the  purposes  of  this  research  eight  BLOCKS  per  ocean 
period  were  considered  adequate  coverage.   Each  of  the  eight 
BLOCKS  was  said  to  be  at  a  different  PHASE  lettered  A  through 
H.   Every  eighth  BLOCK  would  be  at  the  same  phase.   Again  if 
the  ocean  were  monochromatic  and  the  spectral  density  peaks 
in  each  PHASE  averaged,  one  average,  corresponding  to  the 
strongest  part  of  the  signal,  would  predominate  and  the  other 
PHASES  of  lesser  amplitude  would  be  distributed  sinusoidally 
around  it. 

PERK  1A  and  PERK  IB  are  highly  versatile  programs  with  a 
number  of  options  available  as  summarized  in  the  following 
table  (Table  IV) . 
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TABLE  IV 
SUMMARY  OP  PERK  1A  AND  PERK  IB  PARAMETER  OPTIONS 


Total  number  of  A/D  data  points 
Number  of  data  points/BLOCK 

Delay  between  BLOCKS 

Number  of  BLOCKS/ocean  period 
=  number  of  PHASES 


up  to  8192 

from  16  to  512  in  powers 
of  2  (2N) 

0.01  to  999-9  msec  with 
accuracy  to  0.01  msec. 

2  to  16  in  powers  of  2 


The  program  also  provides  a  unique  number  for  each  run, 
based  on  day,  month  and  run  sequence  of  that  day  (i.e., 
#15101  is  the  first  run  conducted  on  15  October) ,  and  also 
prints  out  input  parameters  such  as  sampling  frequency,  prior 
to  displaying  the  output. 

It  should  be  noted  at  this  point  that  PERK  1A/1B 
were  tested  on  a  known  sinusoid  of  2.5  Hz  and  found 
satisfactory.   In  that  test,  groups  of  8  A/D  samples  separated 
by  0.4  seconds  were  taken  at  80  kHz  and  plotted  versus  time. 
No  visible  discrepancy  existed  over  24  sinusoidal  periods. 

To  recapitulate,  the  programs  PERK  1A  and  PERK  IB  were 
used  to  take  BLOCKS  of  64  A/D  samples,  properly  spaced 
temporally  to  provide  eight  PHASES  (A-H)  per  ocean  period, 
over  sixteen  ocean  periods.  This  equates  to  128  BLOCKS  or 
16  BLOCKS  per  PHASE.  Samples  were  taken  at  four  times  the 
acoustic  frequency,  or  twice  the  Nyquist  rate,  to  diminish 
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the  effects  of  aliasing.   A  pictorial  representation  of 
this  signal  processing  scheme  is  shown  in  Figure  42. 

B.   EFFECT  OF  SIGNAL  STRENGTH  ON  MODULATION  - 
LLOYD  MIRROR  MODELING 

One  of  the  most  obvious  times  when  signal  processing 
enhancement  is  important  is  when  the  received  acoustic 
signal  is  weak  and  may  be  fluctuating  in  and  out  of  the 
background  noise.   In  order  to  provide  a  common  basis  for 
comparing  varying  signal  strengths,  it  was  decided  to  employ 
different  aspects  of  the  Lloyd  mirror  phenomenon  to  observe 
varying  strengths  of  the  output  signal. 

1.   Brief  Review  of  Lloyd  Mirror  Effect 

As  every  school-boy  acoustician  knows,  the  Lloyd 
Mirror  Effect  is  caused  by  constructive  and  destructive 
interference  between  the  direct  and  smooth  surface  reflected 
paths  taken  by  an  acoustic  signal.   For  a  mirror  surface  within 
a  Critical  Range   R  =  — r—   ,  the  effect  is  one  of  pressure 
peaks  and  valleys,  with  the  peaks  falling  off  as  1/r  and  the 
valleys  going  to  zero  (h  =  receiver  depth,  d  =  source  depth, 

X  =  acoustic  wavelength).   Beyond  the  critical  range  the 

p 
pressure  declines  as  1/r  .   The  OAWF  tank  Is  capable  of 

measurements  both  within  and  beyond  the  critical  range  for 

frequencies  in  the  vicinity  of  13  kHz.   Another  means  of 

visualizing  this  effect  is  by  replacing  the  reflected  rays 

by  an  image  source  at  I  as  also  shown  in  Figure  4  3. 

In  the  case  of  a  rough  surface  the  peaks  caused  by 

reinforcement  are  not  as  great  and  the  coherent  scattering 
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term,  given  by  the  first  Dirac  delta  term  in  the  expression 
for  the  scattered  sound  spectrum  (Theory  section),  predominates. 
Further  the  valleys  caused  by  cancellation  are  now  minima 
instead  of  zero  as  the  incoherent  or  modulated  term  in  the 
S(ft)  expression  becomes  important.   This  effect  is  shown  in 
Figure  44  which  contrasts  the  rough  and  smooth  situation. 

The  setup  used  in  the  following  three  trials  was  essen- 
tially the  same  as  in  all  signal  processing  runs.   The  OMNI 
was  the  source  and  the  LC-10  served  as  receiver  (both  at  a 
depth  of  10  cm.),  with  a  3  Fan  OAWF  sea,  set  to  the  standard 
water  level,  modulating  the  transmitted  20  kHz  acoustic 
signal.   A/D  samples  were  taken  with  program  PERK  1A  at  a 
sampling  rate  of  80.0  kHz  and  analysis  was  performed  by 
PERK  IB.   A  setup  diagram  is  shown  in  Figure  45.   The  rack 
on  which  the  LC-10  was  mounted  was  moved  back  and  forth  with 
the  water  surface  quiescent  (fans  off)  until  a  Lloyd  mirror 
maximum  and  adjacent  minimum  were  identified.   A  third 
position  was  then  chosen  arbitrarily  between  the  two  extremes 
and  designated  "Intermediate".   Runs  were  conducted  at  each 
of  these  three  locations  and  a  summary  is  presented  in  Table  V 
below. 

From  the  output  of  PERK  IB  it  was  possible  to  compare  the 
signal  fluctuations  at  the  three  positions.   At  the  Lloyd 
mirror  maximum,  corresponding  to  a  strong  signal  (larger 
constant  amplitude  coherent  and  smaller  variable  incoherent 
contribution  from  the  surface  scattered  sound),  the  maximum 
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TABLE  V 
TRIALS  CONDUCTED  AT  VARTPTT5~,T  TCyiV  MIRROR  POSITIONS 


Run  No.   Lloyd    Samples/     Delay   No.  of   Total  Ocean 

Mirror   BLOCK  PHASES    Periods  Covered 


10101 

Min 

64 

47.69 
msec 

8 

16 

10102 

Max 

64 

(assumed 
ocean 

8 

16 

10.03 

Inter- 
mediate 

64 

freq. 
2.578  Hz) 

8 

16 

spectral  density  peak  exceeded  the  overall  average  by  about 
2.1  dB  and  the  minimum  peak  fell  about  2.0  below  the  average. 
The  total  range  of  amplitudes  was  about  four  decibels. 
Applying  the  same  comparison  technique  to  the  other  two  runs 
yields  a  swing  of  about  8.5  dB  at  the  "Intermediate"  position 
and  almost  19  dB  at  the  position  of  the  weak  signal  where 
there  is  a  maximum  fluctuation  of  the  amplitude. 

It  is  apparent  that  in  this  model  fluctuations  do  exist 
regardless  of  signal  strength  and  that  they  are  of  greatest 
importance  when  the  signal  is  weak.   This  further  points  up 
the  importance  of  adapting  knowledge  of  the  sound  modulation 
to  enhancing  a  signal  processing  method  as  a  means  of 
acquiring  and/or  retaining  low  strength  signals  in  a  noisy 
or  highly  variable  environment. 
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C.   "BEST  PHASE"  METHOD 

As  previously  mentioned-  .If  t-.Vie.  scattered  sound  were 
modulated  at  a  single  fixed  frequency  there  would  be  a  PHASE 
which,  when  once  correctly  chosen,  would  always  yield  the 
strongest  average  signal.   The  modulation  frequency  is  of 
course  not  monochromatic  as  shown  earlier  in  Figure  4l  and 
in  fact  its  spectrum  has  about  a  0.5  Hz  half-power  bandwidth 
(Figure  26). 

This  not  withstanding  it  was  felt  both  instructive  and 
possibly  rewarding  to  compare  the  relative  strengths  of  the 
eight  PHASES.   Thus  an  annex  was  appended  to  PERK  IB  which 
averaged  all  of  the  spectral  density  peaks  in  each  PHASE 
(every  eighth  peak  in  this  case)  and  printed  out  a  PHASE 
average  as  well  as  the  overall  average  of  all  peaks.   These 
values  for  8  runs  are  shown  in  Figures  46(a-h),  which  depict 
PHASE  average  spectral  density  plotted  versus  PHASE  with  the 
dashed  line  indicating  the  overall  average.   All  of  the  runs 
are  seen  to,  in  fact,  follow  a  rough  sinusoidal  pattern 
particularly  Figures  46(c)  and  46(g)  but  others  such  as 
Figure  46(a)  appear  to  be  at  double  the  ocean  frequency. 
The  Figures  46(a-h)  show  spectral  density  on  a  linear  scale. 

If  the  "BEST  PHASE",  or  highest  PHASE  average  is  compared 
to  the  overall  average,  the  gain  is  In  the  range  of  0.3  dB 
(for  the  Lloyd  mirror  max-strong  signal  case)  to  0.5  dB 
(for  the  weak  signal),  which  is  hardly  anything  to  write 
home  about. 
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In  summary  the  "BEST  PHASE"  method  applied  innocently 
over  16  ocean  periods  is  seen  to  offer  some  small  advantage. 
However,  because  of  the  broad  band  modulation  of  the  scattered 
sound  modulation  and  the  drift  of  the  best  phase  over  many 
ocean  periods  the  gain  realized  is  minimal. 

D.   "MONDAY  MORNING  SIGNAL  PROCESSING"  METHOD 

The  next  approach  began  with  hand  plotting  the  spectral 
density  versus  time  information  presented  by  PERK  IB  and 
examining  it  to  see  if  any  combinations  would  become  apparent. 
Samples  of  the  plots  are  included  in  Figures  4l,  47  and  48 
and  show  the  fluctuations  of  spectral  density  with  time  for 
the  Lloyd  mirror  minimum,  "intermediate"  and  maximum, 
respectively. 

To  analyze  this  displayed  information  a  clear  plastic 
mask  was  constructed  with  a  sinusoid  of  the  same  frequency 
as  the  assumed  peak  OAWF  ocean  frequency  printed  on  it.   The 
mask  was  then  moved  along  the  time  axis  of  the  spectral 
density  versus  time  plot  to  determine  how  many  peaks  or  near 
peaks  corresponded  to  the  ideal  sinusoid  peaks.   Typically 
groups  of  at  least  three  but  up  to  five  would  line  up  prior 
to  an  irregularity  or  "glitch"  which  would  require  shifting 
phase  to  a  subsequent  group.   This  was  not  unexpected  as  the 
auto-  and  cross-correlation  functions  discussed  earlier  in 
this  work  indicated  good  correlation  for  at  least  two  or 
three  ocean  periods. 
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The  above  procedure  was  felt  to  simulate  a  signal 
processor  which  would  identify  a  strong  signal  and  then 
sample  later,  only  at  intervals  which  are  multiples  of  the 
ocean  period,  three  to  five  times  to  gain  signal  enhancement 
over  continuous  integration.   The  results  of  this  signal 
processing  in  hindsight  —  hence  the  Monday  Morning  label  — 
are  interesting  and  are  illustrated  below: 

Referring  to  Figure  4l,  the  Lloyd  mirror  weak  signal 
(large  fluctuations)  case;  if  the  BLOCKS  indicated  by  arrows 
are  chosen  and  averaged  as  follows: 

Group  BLOCKS  M2SP  Ave. 

1  5,13,21  ''44. 70 

2  39,47,55,63,71      25.20 

3  98,106,114,122      42.55 

These  groups  of  BLOCKS  correspond  to  identifying  initial 
strong  signal  (i.e.  BLOCKS  5,39  or  98),  then  sampling  at  the 
ocean  frequency  for  only  three  to  five  peaks  further.   The 
payoff  becomes  apparent  when  these  Monday  Morning  averages 
are  compared  to  the  overall  average  of  19.6  which  yields 
gains  of  +3.6,  +1.1  and  +3.4  dB  for  groups  1,  2,  and  3, 
respectively. 

Gains  in  the  case  of  the  modeling  strong  signal  (weak 
fluctuations)  ranged  from  +1.0  to  +1.4  dB  and  those  for  the 
Lloyd  Mirror  "Intermediate"  response  were  in  the  +0.4  to 
+1.1  dB  range.   The  interested  reader  Is  invited  to  verily 
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these  values,  or  try  to  improve  on  them  with  other 
combinations,  by  marking  a  ruler  or  like  object  in  intervals 
of  0.4  seconds  and  sliding  it  along  Figures  4l,  4?  and  48. 
The  advantage  of  Monday  Morning  Signal  Processing,  is 
obvious  in  the  weak  signal  environment  (where  it  is  most 
needed,  of  course)  and  is  at  least  worthwhile,  though  not 
vital,  when  observing  a  stronger  signal.   The  author  will 
consider  modifications  to  this  technique  and  suggestions  for 
real  time  analysis  in  the  final  section. 

E.   SUCCESSIVE  AVERAGING 

In  an  attempt  to  validate  statistically  the  "seat  of 
the  pants"  notion  that  three  to  five  ocean  periods  is  the 
maximum  extent  to  which  sampling  at  multiples  of  the  ocean 
period  can  be  profitably  employed,  a  system  of  averaging  was 
developed  and  appended  to  PERK  IB  (which  became  PERK  1C). 

The  spectral  density  peaks  were  first  divided  into  eight 
PHASES  (A  through  H)  of  sixteen  peaks  each.   These  PHASE 
groups  were  then  blindly  averaged  in  ensemble  group  as  follows 

Taking  PHASE  B  for  example, 

.   .   ,      B1+B2    B2+B3          B16+B1 
Average  by  two's    — - —  ,  — 3— ^  ,  .  .  .  ,  5 


.   ,,     ,    B1+B2+B3    B2+B3+B4          B16+B1+B2 
Average  by  three's   = ,  ^ ,  ...  ,  5 


etc . 
Up  to  Average  by  eights 
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These  averages  were  then  compared  with  the  overall 
average  and  the  number  that  exceeded  It  noted.   Rather  than 
confirming  the  initial  hypothesis,  however,  -the  results  were 
contradictory;  indicating  that  groups  of  three  to  five  were 
most  profitable  in  some  PHASES  while  averages  of  five  to 
eight  were  most  beneficial  in  others.   In  fact  one  PHASE 
indicated  that  it  made  little  difference  how  many  were 
averaged  —  that  all  combinations  were  equally  profitable. 

The  author  remains  convinced  however,  that,  though  it 
has  not  been  statistically  proven,  when  compared  with  blind 
averaging,  the  Monday  Morning  technique  demonstrates  that 
averaging  more  than  three  to  five  peaks  in  a  group  is 
counterproductive . 

In  summary  the  blind  Best  Phase  method  provides  only  a 
minimal  gain  of  about  half  a  decibel  when  applied  to  a  weak 
signal  while  the  Mondary  Morning  technique  provides  nearly 
a  +3  dB  enhancement  of  the  signal  processing  for  very  weak 
signals  if  sampling  is  limited  to  three  to  five  ocean  cycles 
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VI .   FUTURE  RESEARCH,  CRYSTAL  BALL  GAZING  AND  OTHER  MUSINGS 

The  author  feels  that  the  research  described  in  the 
foregoing  pages  has  only  scratched  the  surface  of  a  viable 
means  of  acoustic  signal  enhancement.   Some  basic  parameters 
and  features  of  the  modulation  of  surface  scattered  sound 
have  been  identified  or  confirmed  and  in  some  cases  loosely 
quantized.   These  ideas  may  be  of  immediate  use  to  signal 
processors,  however  additional  research  is  felt  to  be 
indicated  in  the  following  two  areas: 

A.   Further  investigation  of  the  phenomenon  for  additional 
clues  which  might  increase  its  value/utility  and  prove  more 
rigorously  some  notions  presented  herein. 

The  initial  concept  of  identifying  a  signal  peak  (or  loud 
point),  then  sampling  later  at  multiples  of  the  ocean  period, 
which  at  this  point  might  be  charitably  described  as 
simplistic,  could  be  more  effective  in  a  modified  form.   If 
a  peak  were  identified  and  the  next  sampling  interval  were 
sought  within  the  interval  T±AT,  where  T  is  the  ocean  period 
and  AT  is  related  in  some  unknown  way  to  the  bandwidth  of 
the  ocean  spectrum  and/or  the  rate  of  "attenuation"  of  its 
autocorrelation  function,  greater  improvement  might  be 
realized.   Of  course  the  interval  of  sampling  would  increase 
by  additional  multiples  of  AT  if  the  starting  point  were  not 
reset  at  a  peak  after  each  group  of  samples  was  taken  and 
the  peak  identified.   (See  Figure  '(9.)   It  is  felt  that  the 
previously  mentioned  limit  of  three  to  five  sampling  intervals 
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Figure  ^9.   Improving  on  the  "Best  Thase"  Processinr  Method 


98 


before  the  diminishing  returns  point  is  reached  may  also 
apply  here. 

It  might  be  even  more  profitable  to  take  samples  at  a 
rate  of  several  BLOCKS  per  ocean  period  over  some  observation 
interval  (greater  than  the  ocean  period)  as  done  in  this 
research  and  then  to  compute,  sort  and  rank  the  spectral 
density  peaks  by  magnitude.   Once  so  arranged  the  following 
options  would  be  available: 

1.  use  only  the  highest  peak  in  each  observation  interval 

2.  use  only  the  highest  peak  in  each  ocean  period  and 
average  over  the  observation  interval 

3.  "use  the  five  (or  any  other  arbitrary  number)  highest 
peaks  observed  and  average  as  in  2 . ) . 

4.  use  only  those  peaks  which  exceed  the  overall  average, 
which  would  also  have  to  be  computed,  by  3  dB  (  or 
employ  some  lesser  standard,  such  as  2  or  3,  if  none 
did). 

As  an  aside  it  should  be  noted  that  when  speaking  of  the 
"highest  spectral  peak"  it  is  implicit  that  the  highest  for 
each  frequency  observed  is  computed,  compared  and  utilized. 
This  research  was  conducted  with  a  single  acoustic  frequency 
but  future  trials  might  include  several  tranducers  operating 
at  various  frequencies  and  signal  output  levels.  A  BLOCK  of 
more  than  6H  samples  may  be  required  in  this  case  to  provide 
the  desired  frequency  resolution. 

Another  possibility  might  be  to  continuously  process  for 
most  of  the  "loud"  half  of  the  ocean  period  (5  seconds  or 
so  in  the  actual  ocean)  ,  then  process  and  use  the  received 
frequency  spectral  information  to  image  correlate  during 
the  quieter  half  cycle.   Again  the  three  to  five  period 
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profitable  limit  is  felt  to  pertain.   This  method  might 
provide  even  better  results  ^^-'v/"t:,~.:>~.  obtained  by  active 
sonar  image  correlating  as  the  changing  doppler  effects  and 
attendant  correlation  mismatches  would  not  be  a  factor  over 
such  short  time  spans. 

B.   The  second  major  area  of  suggested  study  would  be  to 
spectrally  analyze,  plot,  and  employ  the  Monday  Morning  Signal 
Processing  Technique  on  actual  passively-recorded  submarine 
tapes  to  see  if  gains  comparable  to  this  research  are  realized, 
Additionally  the  modulation  frequency  spectrum  and  autocorre- 
lation functions  could  be  determined  and  compared  with  both 
the  locally  observed  wind/sea  data  and  the  results  of  this 
work. 

In  summary,  there  is  no  doubt  as  to  the  fact  that  sound 
scattered  from  the  ocean  surface  is  amplitude  modulated  at 
the  basic  surface  wave  frequencies.   The  ability  to  estimate 
the  ocean  wave  frequency  and  spectrum  is  of  course  vital  to 
taking  advantage  of  this  phenomenon  and  has  not  been  addressed 
to  this  point.   The  author  feels  that  ballpark  estimates 
could  be  arrived  at  on  station  through  a  knowledge  of  fetch 
and  surface  wind  speed  and  duration  and  assuming  a  Pierson- 
Moskowitz  wave  height  spectrum.   Further  it  is  felt  that 
meteorological  organizations  such  as  Fleet  Numerical  Weather 
Central  could  supply  this  information  if  tasked  to  do  so, 
so  that  signal  processors  would  be  provided  with  up-to-date, 
tactically-useful  data. 


I  0(1 


Though  other  mechanisms  such  as  internal  waves  and  tem- 
perature microstructure  also  contribute  to  fluctuations,  the 
surface  scattering  is  of  particular  importance  for  near- 
surface  sources  and/or  receivers.   What  remains  in  the  author's 
opinion  is  to  find  the  means  to  most  effectively  use  the 
various  ocean  modulation  effects  to  enhance  signal  processing. 
It  is  felt  that  clever  signal  processors,  if  properly  in- 
formed, could  devise  numerous  techniques  far  more  subtle  and 
comprehensive  than  the  rather  simplistic  methods  suggested 
in  this  and  the  foregoing  sections  and  that  even  greater 
gains  might  be  realized.   It  should  be  emphasized  that  use 
of  these  techniques  is  not  contended  to  be  the  miracle  cure 
for  the  Anti-Submarine  Warfare  problem.   The  author  does  feel 
however  that  ocean  modulation  phenomena,  if  used  to  best 
advantage,  will  provide  part  of  the  answer  required  to  cope 
with  quieter  next  generation  targets. 
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QUI  •  1 1  m  >    =   rpy<  l .  n ■■■) 


in  ■; 


4?l'l 


out  <  ? -  v  *  >    =   apy  <:?:  •  k"IO 

OUT  .  l  iV'K- 1  >    =    RPY  ■  1  •  KK- 1  > 

0'  IT  <?. .  M  -  1  •■     =    APY  <£ .  K-  -  r> 

K    =    f^ 

Tjn    4;?0    1  =  1  jHV.P 

fc'K    =  *  /? 

TEMPI     =    FLOPT-rifiPYtf  ,K'0 

TEMP?    =    FLOhT  c-I  APY  <-  5  >  fc - 1  ":■  "> 

nijTr  1  ,  kf  i     =    P,PY> '«  .  k>:  > 

OUT  <•-; .  i'i-  *     =    ftPV  ■  '£• .  V  IC> 

PPY<1.*">     =    TEMPI 

ftp 7  ip  ,  |t';,      =      n  _  0 

pc-v  .: 1 . »-  -  ]')  =  TEMP£ 
PPY  <"£«  K-i ')  =  0.  0 
V    -   *:-£' 
CONTINUE 


4Qfl 


IF  aL-n  7*6-7 


HO  3  I=1«NV.S 
IT  =  I +? 

up  i tf  ( ■'< .  a n£)  c apy  < i -  mi  .  m=  t  « 1 1 : 

POFMAT     i'SF7.  IT' 
CONTINUE 

'.'PITT    '"zusro 

PFhTi     'S.o^O)      ICO 

IF     CICO.E'O.  n.)    i?D    TO   8  0  0 


A5  fi 

t'-1  I 


IX    =    -1 

CALL    FOUP 1 < flP Y . N V « I X ) 
I.IPITF      f?,  87  m 
PEATj    (5«35<V>  ICO 

IP    CICO.EQ.  i"1''    60   TO   A  on 
IF     ('IP.  Ffi.  (1  .     hD    tc    6.1 1 
i.ip  i  tf     •  ■'■■  ■  r.  1  0 > 

prjIPMPT     'lif'»'.  5HPSP"  15"-  •  4HFPFO,  t.X«5HRNGLE/"' 
p<;   =    n.  u 


output  i  1ST 


1  ,-■  no.  COMPLEX  point*  +  s 


s  n  1 

50? 
5  fl  i ' 


SO' 


pn    .:•     J  =  j.  IN',,',  ip-jfl 

I P  C  AP  Y  '  !  •  i  •  ,  ME .  0  .  0  '  i"  0    TO    5  r  0 

ANSI  r     =;;•;■  p.  fi 
■"-.n    i  n   f,  •  ij 

fcNAI    F='l.    fl 

GO   to    53 f i 
ANGLE=?n. 0 

^0    TO    *"•?''' 

MtluL  F=ATAN  '  PP-'.'  '?.  I  •>    'APY  <  1  <  I  ">  > 

]C      ,HtY  .  1   .   |  ,  ,  ,,t  ,  ,,.  ,;,  ,      (ZQ     JQ    S<15 
AMRLE     =     'Mm. I   |    *  111.  n 

'"0  in  f  "  n 

I  it   /rti.",' .?  .  i  •  .  i.r- .  n .  !-•-,    r.n   ip   5  ->,-, 

PI T"    -    APY  ■  «  •  I  ■  ♦  ■  :-     <  i  .[...•;.■,' '.'-■.  r  •  ■<■<'<-'■'.' 

II  <   lh   .(".I'  >      AO      F'l     61  < 


I'' 


10 'I 


WRITE     '3'PQO'     PSP«P" 'ANGLE 
61?         ftPV*  t «  I  •     =    P~T> 

fipv<?,  n   =  flna  e 

tirtii         FOFMFiT     .-?.v.  Flf  .1  .  r:.F14.3«3X»FS.3> 
IP  <PS.LT.TE>     GOTO    £0 
RS=fl .  I'l 
GOTO    ?  ...'"' 

2  COMTTNUP 


r 

r 


IIP  1 7E     '3-  P7fl  > 

IP     fTCQ.EO.  i"''    <5D    TO   S00 
lc.  CONTINUE 


r 
f 


r 


PS  =     0 . 0 
l.icjTF    (-i,A*.Qi 
4^.0         PDPMfiT     .vr?v-4HFPFr>,?X,  1  h*\  OLOG  'P7  Dl-PST'c  '>  -3X-  13HP.NGLE1-RNGLES  -O 

T'O   430    I  =  t  >  IHV,  IPSP 

I  F     fflPV  <  1  '  I  >  .  '€  .  0  .  TO     GO    TCI    4S  0 

P5H  =    99'^q  .  n 

GO   TO    4"M 

45 0  P "  P  =    1  0 .  0*FC 0G1  0  'OUT  •  1  .  I  >  /flPY  >  1  .  I  >  > 

451  hHi"LF    =    niJT  •£'  '  P-fiPV  ■  ■"£•  T "• 
MCITF     ''3«46D    PS.PSP-.fiNGLE 
IF     -rPS.l.T.  TB?     GO    TO    440 

p.:    =   0.  0 
GO   TO   4?0 

44  f.  ?-:     -     PS+PES 

4?n  CONTINUE 

46 1         FORMRT    < 3V >  F 1 4 . 3.  3X • F 1 6. 5? 3X «  F 1 3 . 3"> 


GO  TO   3  00 
STOP 
EH? 


ID'.. 


I. 

r  ppc-^lfi      ?_t    RUGU'T    l'='74  rcH 

0  pppRPOM   TQ   SPMPLE'POIJSEjPNB   MPITE   PRTfl   dh   cassette 

r  P    SPECIFIEP   Nl.WPEP    OF    POINTS    RPE    TPtEM   FPOM    flPC-r    IN    P 

••'  SPECIFIED  T-i.  nr>     :ize    "."TTH    h    SPECIFIEP    PELRY    BETWEEN 

r  FPr.H  ftn> 

r  INPUT    pPTP: 

r:  PIJH    =    TPENTIFYIUG    HIJMFEP     flMTEGEP^ 

C  N   =    PT3I.IEP    OF    ?    EOUPl    TO    PPPRY    SIZE 

C  MS    =    FftuFF    OF    .?    EC'UPL    TD    BLOCK    SIZE 

f  SF    =   ShMPLF    FPFnijFHCY 

f:  HP    =  Hl.MFFc'    OF    FWfi"~ES 

C  T'FLPV  =  T'Fl  hv  pCTWFFM  FLOCKS  IN  MILLISECONDS  ♦  100 

.-  1  =  YES.  0  =  NO 

C 

I  NTFGE  P*p    N .  TNP .  NS  «  PPS  '  PEL  PY «  CON 

INTEGER*?    I  ,  INPX-  IPPYjKENTJj  1 1  •  M. 

INTEGER*?.    HP<    I.  JJ 

INTEGER*?    htote 

INTEGER*?  PIJH 

PIMENi  ION    IflPY  e8£001 

KENH    =    ?.£7G7 
f 

f:  REQUEST    INPUT    PPRPMETERS 

r 
inn         i.ipjtE    .?.9  0rt,> 
QOfi         FGPKPT     t'SOHN    =  t     fl£)     <00=ENIO     > 

pfhP    ''f,.?iM>     N 

QflJ  C-flpMpT       ,'If, 

IF     i.'N '■•      1  OO'c  Ol'l'  1  1  0 
110  IF     (N.CT.  I?)     iT    TO    100 

TMF    =     ?**N 

r 

ISO         I.IPITF     <3'  ■?*)?*> 

■up?        pppmht    <l:fWNS    =    "'    'IS')    > 

ppH-n    <"5,9no    NS 

IF     nt~ :•'■<     1  £fU  130»  1  30 
130         IF    fNS.GT.9"i    GO    Tn    iso 

TF    nr,iT.4'    GO    TO    ISO 

ppt;    =    ,-.♦,•.; 

MTOTP     =     TNP/'PP?  j 

IF     tMTGTB.  GT.T-l.S'i     GO    TD    1  Or. 
TNP    =     NTGTB*PPS 

l"  I 

I.IPJTF     •  ] ■■.  --rr.\ 
■-tp':  FPPMPT     <lf.WpFI.OY    =     '■      <  \f,')     > 

PF.fip     •t.'iii4.     DELRY 
904  FGPMP.1      .IS,  ^    j 

.   ' 

I.IPITF     <"5'9ns^  ■  .' 

305         FOPMP.T     f£SHR.«N  NO.     =    ?     <I5'<     <.'IT"-1RP>     ■> 
REPP     .-'-, ,  ---.n  a  .     pi.in 

I.IP1  if     ■  5,«rt«) 
90?  FnPNfiT     .I.IH'F     s     v     ,frc:-p     ', 

pppn     .  c, .  Ql  ft".        (  .' 

•^I  0        FPPMPT    -Ff.  i  >  .-'..'■. 

1  4  0  MR]  11"     •  .-;.  -Ml' 


r 


r 
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r 


PPPP      rS.'Jni'      HP 

jf    <■  CNTOTE --MPWNP.^f.ntoTB''    GG   TD    I4i1 

WRITE    f?«?0*:'1 

9(1*  FPFMfiT     .'1  1HCDNTTNI.C     '""    1 

PFPTi     •. "? .  ■=•  fi 7" :.     CQM 
907        fopmht    t  1 1  '■ 

IF    i'CON.FD.  mi    RQ    TO    IOC' 


I  =  TNP'PPS 
TMP  =  I*PPS 

CLEfiP  THE  P.PPRV 

rio  i5i"'  1=1 » 9? no 
ippv  en   =  t'END 

150        CONTINUE 

r- 

i"        GET  I'RTP  FPDM  PUC  CHANNEL  C 
r 

CPLL    SfHMP  •'  I RPV .  TMP •  FP S .  IE LfTO 
r 

r  OMTPUT    IihTP 

f 

I.IPTTE     C?«  ^S"**     PUN «TNP,PPS« DELAY,  SF.NPjNTOTF 
??  i":         cC.PriPT     'M  T  £ .  FS .  1 ,  2.1  f.-"> 
90S         FCPMhT    (SI  61 

tin    l^n    I  =  1.TNP.8 

II  =    1+7 

I.IP  I TF     .:? .  9  h:?>     <'  I  PRY«'  M}  .  M=  I «  I  r- 
1f-M"i  CONTINUE 

MPITE     .a,9  0fi>     KENP 

■?o  to  ion 


i 


?00         i.iPITF     •S,.9n9i    KENT) 
STOP 

i  EMT' 
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VC--.  "   •'  ''•":■'  '   '.     .'•'. 

:    C  ■'■•■•'  PEPTIC      15    OCTOPEF     W74  OCH 

,    C  PROGRAM  TO    RMRLYIE  PLDCK    DRTR 

C  • 

C  •  ALL  INPUT  PAPAMFTEP'  AHP  PhTR  APE  DM  R 

C  •'  ■•■•'  CASSETTE  GENERATED  PY  PFPK1R 

c 

DIMENSION   RPY'fP.^.12')  •  PHA'E  <?.,  51  3^  .RPSIi<£.  16) 
f: 
C  GET     I MFI.IT    FRPRMETEPS    FPtlM    TRPF 

r:  : 

100         REAP     (?)?fin^     NPIJN.NTNP.NPPS.NDELRY.  SF.MP.NTOTB 
900        format    (4I6.fr. i.si£>     . 

IF     'HPI.IN.EO.  3?7ir-7">     50  TO    4fin 
UPITF     (3.902* 

902  FOPMAT     IT9WTPIHAN    ?     ) 
PERU     >:5.9n3*     IT 

903  FOPMAT     'ID 

!     •  WRITE     (3»  930* 

!        930         FORMAT     fp^HFWri    ?    > 

J  PEAD     •:5-?n?N     IE 

IF     ciE.Ff'.l'     GO    TD    400 

I.IPITE     (3.94  0* 

940  FOPMRT     <1*HTEST    PRINTOUT    ?    * 
PEAP     0?.903*  ITPO 

I.IPITE     (3.941* 

941  FOPMRT  (?4H INTERMEDIATE  PRINTOUT  ?  * 
PERU  (5»903*  IPO 

i    c 
;    r 

j  liiPITF    io»931> 

991         FORMAT     '  ' 4 sn ♦-!>»♦♦  ♦♦♦♦♦*♦•♦•♦»♦*♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦*♦♦*«♦♦♦<>♦»♦♦> 

I.IPITF     (3«913*     NRIJN 

913  FOPMRT     "V3HRUN? IX. I5>  . ' 
I.IPITE     (3«914*      "F 

914  FORMAT     (llH^AMPLE    FREP,  1X«  F8.  1> 
I.IPITE     ■3.915''     NpFLRY 

915  FOPMRT     (5HPELRY)  1  V.  I5> 
I.IPITF     (?  j  93 11 
MPITE     (3.932* 

93?         FORMAT     <•* 


r 

r 


DO    270    np= 1. MTOTE 


_ 

PO    110    1=1 . NPPS. 3 

II     =     1+7 

PFAD    '  (?.  ,r,0l  '.     (AF  ',  •  i 

1 1  0 

CONTI  Nl  if 

Mill 

cn^l'WT      . «  F  £ .  II  . 

do  i?o  i=i.tpp: 

RPV  .  ?,  I  <     =    0.0 

130 

CDNT1  HUE 

c 

r 

IF    ''IT.fo.  o*    GO   in 

r 

r 

j= i ■  i  r 


C    =    c  .  II-  Fl  nRT  .  HPPC> 


m,,-; 


*  i-j'V'-r.'-  .  g-.cr   r  ■•.:•■     f\  -,'   ■      &■  .:.    -•-  ".       ■•":  '-  •'-••■ 

::-V^:  «'-  C  ';...  .    ,;.;„■ . 

.'•- ..->&***••  /-DO   140   i  =  i,npps  .,       '.....    

.".«•• ':-:'.'-  -flRY«fl . I  >    =   PPVfj.I}*?  •    • 

.1 ;;V»>. •'/  IF     CR.GE.  l.fT<     GO    TD    130  .    •  ■ 

.'"•v.'..-«  •".  B- «   B+C  :.  .■''.''*      .'.'  ■■/.  ";,■'.■':<   ■"■'... 

■■}■■    -'■''  ;..-en  ttj  135      •  -■  .•■  v  -:\f :'     ■'■■■'.   .'."-'J.-.'.-  • 

•'•.  130         B   =    E-C  "' .    ■'   '    ;.":' 

"'  135;"    A   =   fi+C  •      '"  ■",  .  •'  '•' 

*  i40     ■continue  V    ■   -'■■-.;' 

~c  ■     '  -  ■  ...:;•'•.'■    ■    . 

c " 

150  IF     nPO.E<?.0>     sn    TO    410 

I.IPITF    <3j904-> 

904  FOPMAT     f  '1  OX'  ^("•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•> 
I.IPITE     '.'3*905*>     NPI.JN.NP 

905  FDPMPT     O'lOX.SHPIJN,  IV.  I "5 , "I  0X.5HPl.DCKj  IX*  14) 
MPITE    >'3«90<V> 

r 
r 

4)0         PES    =    SFrFLDPT  .TMPPS> 
TB    =    SF/p.O 


r 
C 

C 

r 


C 
C 


PHASE  Tl.NFI    =    0.0 


IX   =    -1 

CALL    FDUP1  CflPY.NPPSjIX'* 


IF    ncO»EQ.0'»    150    TD    4£0 

MPITF     c3.90'O 
90G     .    FORMAT     •■  'lfcX.SHPSDj  13X<4HFPE0«6X>5HPMGLE-''> 
4£0        PS    =   n.n 

indx  =  wpps^+e 


PD    £60     I  =1  .  I  MUX 

I F    < BPV a  .  1  '>  . HE .0.0  >    GD   TD   £  1  0 

I P     <APY  (£  iD)     1 8  0  •  1 9  0 '  £  0  0 
ISO         ANGLE  =    £70.  0 

GO    TD    £3  0 
190         ANGLE    =    0.  0 

GO    TD    £3  0 
£0  0         RNGL  E    =   on.  n 

GD    ID    £30 
£10  ANGLE    =    HTPM  •  OPV  i'£.  I")   -fifty  ^l.  I'm 

RNGLE    =    Ht.i-Lt  ♦57.  £*5~"7,:i51 

IF     .hPY.'I  «  I  1  .  i"F.  0.  0>    GC    TO    ££0 

RNGl E    =    RNGl F  +  IS'O.  0 

RP    "m    £3fl 
?£fi         ]F     fpPY<£.  T>  .  GE.  0.  01    GD    Tn    £30 

HtiR  E    =    RNGLE+3^0. 0 
£30         PSr    =    PPY<  1  .  I  ■  ♦  RRY  •  1  .  1  ■  +hPY  ■  Lr',  ]-<  ♦Htv  •  ■£.  D 

IF     iPSP.LT.  FHh"F  '  1  .  MpO     GO    TD    £40 

PHm-E  .  1  ,  tip  ■     =    P"P 

PHASE  1  :.  <  tip-.    =   p? 

=■4  0  IF      'IPO.Fn.iV'     |?D    TD    4?0 

MR]  7F     •  : ■  •  '•nr '•     F'Tu  F  ■  .  RNGL  E 


11)0 


r-^-r-^ ., 

'.90?'    .'FOPMhT     C?X.Ft6.  1i"-iX,F14.3,3X.FR.3> 
•  :."'43tl         IF     'Rv.LT.TP)     GD    TO    350 

';•",;•    P3  =  o.o 

';;  ''-■.•.*•'■       GD   TD    260 
:  5r50        PS    =    P3+PES 
260    '  'CONTINUE  '  : 

'■-'.    WRITE     C3.932">  ■■    •.   . 

WPITE    <3»932> 


V  _ 


|?70         CONTINUE 

W»ITF     <3,942"> 

942  FDPMhT  CSV.f.HBLDCK.  1£X»7HMRX  PSDj  13Xj  4MFPEC» 
TO  44  0  1=1  . NTOTB 

WRITE  <:3.943>  T  .PHH'Ea  .  P  .PHASE  <2i  P 
44  0    CONTINUE 

943  FDPMhT  f 4X . 1 3 . 3X . F 1 6 . 1 . 3X . F 1 4 . 3  > 
C 

r 

PXV   =    FLCPT'MTOTP-HPI 

no   31 0    1=1 « HP 

psn  =  o.o 
ii  =  i 

FPEO  =  PHASER-  P 
HD  23  0  J= 1 , NTOTB . MP 
PSIi  =  PHASE ' 1  -  1 1 ^ +PSTJ 

if   ''iTPo.EO.r^   gd  to  4^0 

I.IR  I TF     <T3  .  944  >     PSIi  •  PHP  3  E  <  1  .  ID 
Q44  FDPM9T     (SFlb.P 

450  IF     I'FPEO.HE.  PHASE '2.  11^     GD    TD    300 

II  =  II+NP 
380    CONTINUE 

RPSDcit  p   =  pstvihv 

IF  -riTPD.EQ.n-.  GO  TD  460 
MPITE  (3 .944')  Hp~- p.  i ,  i) . DIV 

4A0    II  =  I 

SUM  =  0. 0 

IiD    29  0    .1=1  . NTDTP.NP 

PSH    =    PHASE  ( 1 «  II?'  -RPSTJ  •:  1  -  I ') 

sum  =  ji.iM+  >rpsri*psn> 

II  =  II+MP 
290    CONTINUE 

APSTK2»  P  =  SGRT  <3UM.-"HIV) 

GD  TD    310 
3i"i0  I.IPITF     .'3.ci  ri:r;-,      j 

■ans  FnPMRT     ^22HFPFl7i    DIFFERENCE    PHASE     .  12^ 

:-:i0         CONTINUE 


I 


P3P  =  n, o 

WRITE  <3«9'<D 

i.ic  T  TF  .  ?.91  0-) 
910    FORMAT  •   :•:.  chf  hH'F.  1 ,: : '«  ri-M'.'rtHKF-  7".  !  2H"  TfiNPARI1  T'EV) 

PO  320  1=1 .MP 

WRITE     <3«909>     I .  PP^ri  •  1  •  P  .APSIi-'2.  p 
9fl?  FDPMHT     .-sv,  12.  3V.F  1 6  .  2  .  3>> F  1  6. 2  > 

p:  p    =    hF  :t>-  1  «  P  +P  Sfl 
320         CONTINUE 

MVr,  =    F-T-'FLDhT  <NP^ 

■  i  IM  =     0  .  n 
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.-.-.-•.  Iirj   :-r?n    i=i,mp  -~-  •  .■  '•"."; 

y  '_-■';  PSD   =   fiP^TKl.n-ftVG  .; 

r'    ••  '»•  S1.IM   =   tLH'W- •.'P<TP*PrD"> 

•  330  CONTINUE 

;'-''-  SB  =  'SOPT<'Si.»1/'FLOP.T''NP>>  '    "  .'   .. 

■«:'-%.  '-"■  WP1TE  (3,931> 

"7""  WRITE  •'?.9m  fl'-n; 

i-911  FOPMRT     </3X,  lPHnVEF^LL    PVEPPGE    =    >F1<&.£> 

'•  WPITE     <3'<?\2'>     SP 

:"  91?  FDPMflT     <3X,£9HDVEPflLL    STPNTtRPP    PEVIPTION    =    .F16.2> 
WPITE     ^3j931> 

c  • 

HBP    =    MTtTTB-NP 

mpite    <3«93£> 
l-IPITE     •:  3  -  *  ?  1  > 
Tin    "i?fi     JJ=1  .  NP 
MPITE     -'J.siSrn     J.l 

950  FDPMflT     O-IM.^HFHPAF.  IX.  I£> 
I.IPITE     <'3.955"< 

955         FDPKhT    f  .•*r.y.=.wCYCLE.  17X>  ?HPSP> 
NH    -    JJ 

pn  ^05  n=i.n?.p 

MP  I TF     ■:  ?• «  95£>     N .  PHP"  E  <  1  «  NN> 

MN    =    NM  +  t>P 
I     956         FDPMRT     '■?<-  12.  4X«  Fl  £.  £> 
'     505         CONTINUE 

WPITE     ^3.957"> 
j"  957         FOPMFiT     c"  -'I  X .  5HSTflPT'> 

WPITE    i':"::.951") 

951  FDPMPT    ■nx«5HCYCLE«  1  ijx.  4HhVR£.  1  0X.4HPVP3'.  1  0'*:,  4HPVG4.  1  0' -',  4HhVi?5> 
VPITE    <3.?5£':> 

<=i5c  FDPMflfT     '  3 OX .  4HRVhK «  1  OX .  4HfiVG7 «  1  OX .  4H0VG3) 

I  =    JJ 
K   =    1 

'  PP   510    J=t.NBP 

"I  I    =    I +NP 
IF    d  1 .  6T.  NTOTP.}     II=.JJ 

suns  =  phhs  e  <  l .  i  ■>  +phpse  ( i  >  1 1  "> 

ii  =   ii+np 

i^   ni.RT.  N-rn-m    ii=jj 

SUM3    =    SI.W£+PHP?E''1«  II) 

II  =    II+NP 

IF    fll.GT.NTDT?'*     II=JJ 

-|.;m4    =    -.iim-:  +  phh:f  ■  1  .  II) 

11=     II+NP 

IF     •!  I  .'"T.t<TOT?>     11  =  J.I 

-V.HM5    =    "i  inj+PWHrF  .  1  .  i  i  . 

II    =    II+NP 
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